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Теоретически показано, что мягкая мода в окрестности структурного перехода исходная структура — предпереходное состояние 
— конечная структура, вынужденного как изменением температуры, так и внешней силы, переторможена и, следовательно, ди- 
намика параметра порядка имеет релаксационный характер. 


1. Введение 

Экспериментально обнаружено, что кристалл, 
неустойчивый относительно структурного перехода 
исходная - конечная структура, вызванного как из- 
менением температуры, так и внешней силы, в ок- 
рестности структурного перехода находится в пред- 
переходном состоянии [1, 2]. Для теоретического 
описания структурного перехода и предпереходно- 
го состояния сформулирована микроскопическая 
модель [3], описывающая структурнонеустойчивый 


кристалл как квантовую систему псевдоспинов. 
Причина возникновения предпереходного состоя- 
ния в кристалле в том, что изменение внешнего 
воздействия (температуры, механической силы), 
стимулирующее структурный переход исходная - 
конечная структура, уменьшает площадь горба, раз- 
деляющего минимумы двухямного кристалличес- 
кого потенциала атома. Это приводит к существен- 
ному увеличению квантового туннелирования ато- 
ма и уменьшению асимметрии двухямного кристал- 
лического потенциала. Возникает неустойчивость 
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состояния исходной кристаллической решетки с 
асимметричным двухямным потенциалом относи- 
тельно возникновения предпереходного состояния 
решетки с симметричным двухямным потенциа- 
лом. В результате на верхней границе неустойчи- 
вости (при температуре Т* ) возбуждается мягкая 
коллективная мода - псевдоспиновая волна, пере- 
водящая кристалл из исходной структуры в предпе- 
реходное состояние. При дальнейшем уменьшении 
внешнего воздействия на нижней границе неустой- 
чивости (при температуре Т ) происходит конден- 
сация мягкой псевдоспиновой волны, которая пе- 
реводит кристалл из предпереходного состояния в 
конечную структуру. В рамках приближения моле- 
кулярного поля (ПМП) [3] можно изучить только 
термодинамические (стационарные) свойства 
структурного перехода и предпереходного состоя- 
ния. Поэтому для изучения динамических свойств 
структурнонеустойчивого кристалла, например, 
частоты мягкой псевдоспиновой волны, необходи- 
мо выйти за рамки ПМП и использовать приближе- 
ние хаотических фаз (ПХФ) [4]. 

В настоящей работе с помощью ПХФ показано, 
что при почти непрерывном переходе исходная 
структура - предпереходное состояние - конечная 
структура собственная коллективная мода системы 
псевдоспинов смягчается и спонтанно возникает 
разница заполнения левой и правой ям потенциаль- 
ного рельефа, то есть происходит Бозе-конденсация 
псевдоспиновой волны. Если не учитывать релакса- 
цию системы псевдоспинов в окрестности структур- 
ного перехода, то это приводит к динамике с неза- 
торможенной мягкой модой. Однако, хорошо изве- 
стно, что вблизи структурного перехода динамика 
параметра порядка является колебательно-релакса- 
ционной и описывается уравнением типа Ландау- 
Халатникова [5], причем мягкая мода, как правило, 
переторможена [6]. В этой связи возникает задача 
учета релаксации системы псевдоспинов в окрест- 
ности структурного перехода исходная структура - 
предпереходное состояние - конечная структура. 

2. Гамильтониан системы и основные уравнения 

Поскольку нас интересует отклик системы псев- 
доспинов на компоненту механического поля нап- 
ряжений 0,(0, изменяющую асимметрию двухямно- 
го потенциального рельефа и зависящую от времени 
и пространственных координат, запишем гамильто- 
ниан квантовой системы псевдоспинов в виде [3] 

Н = -X НВД + [(1 / 2)1>Д г + 

І І 

+(1/з)Хй/, ; „д;^+ш,о)]^}. (і) 

І,т 

Здесь 8 х , 8у, 8] - операторы Паули для спина Уг, 
Нщ - расщепление энергий четного и нечетного 
состояний атома в двухямном потенциале, Ы ѵ 
(Щ,„) - константа двухчастичного (трехчастичного) 
взаимодействия псевдоспинов. Вследствие вре- 
менной зависимости Н средние значения спино- 


вых переменных также зависят от времени, поэто- 
му будем обозначать их 8?(1)=8°> і . 

Для исследования динамики системы псевдос- 
пинов, описываемой гамильтонианом (1), исполь- 
зуем стандартную методику построения теории ха- 
отических фаз [4]. Гейзенберговские уравнения 
движения для средних значений спиновых опера- 
торов (в системе единиц, где Ь=\ ) имеют вид 

(*)/* = -і < [|,Я] >, . (2) 

В приближении хаотических фаз, являющемся 
прямым обобщением ПМП на случай задач с вре- 
менной зависимостью уравнения движения (2) 
приобретают вид 

с18 і (I) / сіі = 8;(1)хк;(1), (3) 

где зависящее от времени молекулярное поле к,{і) 
дается выражением 

/г,.(Г) = -[д<Я>, /ВД]. (4) 

Эти уравнения эквивалентны уравнениям для 
свободной прецессии псевдоспина вокруг мгно- 
венного значения молекулярного поля в данном 
узле. Однако в ур. (3, 4) не содержится релаксация 
системы псевдоспинов, которая, как было сказано 
выше, играет принципиальную роль в окрестности 
структурного перехода. Остается, таким образом, 
дополнить (3, 4) членами, описывающими релак- 
сацию. Вначале заметим, что (3) обладает хорошо 
известным интегралом движения 

8 2 = (8 х ) 2 + (8 У ) 2 + (8: ) 2 = (8 0х ) 2 , (5) 

где - термодинамически равновесное значение 
8 х в предпереходном состоянии кристалла. 

Структурный переход и динамика параметра 
порядка в его окрестности - это коллективные яв- 
ления. Отсюда можно предположить, что эффекты 
смягчения собственной коллективной моды и ре- 
лаксации в окрестности перехода также имеют 
коллективный (когерентный) характер и поэтому 
не нарушают закон сохранения (5), как это имеет 
место в теории ферромагнетизма [7]. Следуя и 
дальше отмеченной аналогии, перепишем (3) с уче- 
том релаксации, опуская индекс узла /, в виде 

с18/ск = 8хк + К, (6) 

К = -Л[8 х[8х к]] = -Л[(8к)8- (5°') 2 А]], (7) 

где К - релаксационный член, а Я - феноменоло- 
гический параметр. В [8] было показано, что в низ- 
кочастотном приближении механизмы нелиней- 
ности, с одной стороны, дисперсии и диссипации - 
с другой, могут быть учтены аддитивным образом. 
Используем низкочастотное приближение, выра- 
жающееся в нашем случае неравенством 

со < тіп(® с ,т _1 ). (8) 

Здесь со - характерная частота динамического 
процесса, щ, - частота мягкой коллективной моды, 
г - время коллективной релаксации параметра по- 
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рядка 5'. Принимая во внимание замечание о воз- 
можности аддитивного учета нелинейности, с од- 
ной стороны, дисперсии и диссипации - с другой, 
в приближении (8) линеаризуем (7) по 5’, 5 У , 5 х , С1. 
Введем эффективную статическую восприимчи- 
вость х е/ , определяемую соотношением 

5° = х еГ Ь°, (9) 

где У и /Н - соответствующие равновесные значе- 
ния 5 и /; . Из (3, 4, 9) при 88 /сІІ= 0 следует, что 

Х е/ = -5 0г /©„ = а,- 1 ЩПа> 0 /2к в Т), (10) 

где к в - постоянная Больцмана. Здесь использова- 
но выражение для термодинамически равновесно- 
го значения 5 х в предпереходном состоянии крис- 
талла [3]. Тогда после линеаризации К из (6, 10) на- 
ходим векторное уравнение 

с15/ск = 5хИ-у(5-х е/ Ь), (П) 

у =-Яю 0 5 0 -' »Я® 0 2 /2/. (12) 

Компоненты вектора молекулярного поля, сог- 
ласно (4), имеют вид 

к = (~со 0 ,0,3 0 8 : + І 0 (8 : ) 2 + 0(0), (13) 

где /о— І 0 —8,І^. 

1 ],т 

3. Уравнения динамики системы псевдоспинов 

без учета релаксации 

Поскольку 57=0 (11-13), как в предпереходном 
состоянии, так в исходной и конечной структурах, 
то уравнения движения для компонент вектора 
флуктуации псевдоспина 55, без учета релаксации 
приобретают вид 

ісо55 х = іХК ; 5; 2 Ш + 

3 т 

+&.+2І + О 4 ) 

) т 

ісо55:’ = Н-5ГХ1Л + 2 2ХХ']<55;} - 

3 т 

У т 

-5,°Ж,. - (ХіЛ + 2^ я 5°:}58 х + а,}58 х , (15) 

У т 

ісо58; =-со 0 88 у . (16) 

Теперь (14-16) - это система 31Ѵ нелинейных 
уравнений, связывающих амплитуды флуктуаций 
псевдоспинов относительно молекулярного поля с 
амплитудой зависящей от времени компоненты 
механического поля напряжений. Собственные 
частоты системы псевдоспинов можно вычислить, 
полагая О,=0 и решая однородную систему линей- 
ных уравнений (14-16) относительно со. Представ- 
ленная выше система 3/Ѵ уравнений может быть 
сведена к N системам трех уравнений для заданно- 
го волнового вектора путем преобразования к нор- 
мальным координатам. Фурье-образы средних зна- 


чений операторов флуктуации псевдоспина явля- 
ются искомыми коллективными переменными: 

55- = ^<55 ; ехр(-іЩ), 80,- = ]Го ехр {-ІЩ ). (17) 

Рассмотрим свободную прецессию псевдоспи- 
на в предпереходном состоянии кристалла, достиг- 
нутом за счет изменения температуры, где 5, 0г =0. 
Используя координаты (17), получаем для прецес- 
сионного движения с волновым вектором д систе- 
му нелинейных уравнений Блоха 


ісо88;=Ѵ-881+8СІ-}88 у , 

(18) 

ісо88 у ={со 0 -3-8 0х }881- 


-5 0х 8П--У-851 + дП-\85;, 

(19) 

код 8 х - = - со 0 88 у -. 

(20) 


Однородная система линеаризованных уравне- 
ний (18-20) при йЮ ? = 0 имеет нетривиальное реше- 
ние только в случае, если определитель этой системы 
тождественно равен нулю. Решая секулярное урав- 
нение для собственных частот в духе [4] , получим 

соМ) = 0 , ( 21 ) 

ю 2>3 (?) = ю 0 (ю 0 - ^5°*), 5 0х = (I/ 2)Щсо 0 /2к в Г). (22) 

Решение а>і((7)=0 (21) соответствует продоль- 
ной моде, то есть флуктуации в направлении моле- 
кулярного поля, параллельного 5 х в предпереход- 
ном состоянии кристалла. Два других решения 
й) 2 , з(<?) (22) отвечают поперечным флуктуациям и 
описывают свободную прецессию псевдоспинов 
относительно молекулярного поля. Однако на 
нижней и верхней границе интервала ( Т , Т + ) час- 
тота прецессии становится малой величиной всле- 
дствие ослабления двухчастичного взаимодействия 
псевдоспинов когда 

2© 0 -/^іБ(ет 0 / 2к в Т) — > 0, (23) 

то есть при Т-> Т с . Разумно предположить, что, по 
крайней мере, в пределе Т С -^Т Ш (23), структурный 
переход из предпереходного состояния в исходную 
(конечную) фазу связан с конденсацией мягкой 
псевдоспиновой волны. В окрестности Т с можно 
разложить ш 2 ,з(4о) по степеням Т- Т с 

= (дсо 2 /(сі 0 )/дТ) т=Те {Т-Т с ) = 

= Щ 0 ){Т-Т с )/Т с , (24) 

Щ 0 ) = (со 2 3(д 0 )/4к в Т 2 ) С Ъ~\со 0 /2к в Т с ), (25) 

в то время как ш 2і з 2 (^?о)г= 7:^9) г Д е (7о волновой век- 
тор мягкой (критической) моды. В предпереход- 
ном состоянии при Г,—» Т ш квадрат частоты мягкой 
псевдоспиновой волны, как и в случае мягких фо- 
нонов, является линейной функцией Т->Т С (24, 25) 
в окрестности точки Т с . 

Исходная (конечная) структура определяется 
"замороженными" смещениями в левую (правую) 
яму двухямного потенциального рельефа, отвечаю- 
щими мягкой моде. Параметром порядка, который 
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представляет собой замороженную координату, со- 
ответствующую мягкой нормальной моде, в данном 
случае служит однородная спонтанная разность за- 
полнения (поляризация) левого и правого миниму- 
ма двухямного потенциального рельефа. Положи- 
тельное значение разности заполнения соответ- 
ствует исходной структуре, а отрицательное значе- 
ние -конечной фазе. В предпереходном состоянии 
разность заполнения левого и правого минимума 
двухямного потенциального рельефа равна нулю. 

В состоянии исходной (конечной) структуры 
величина 5 ю *, так же как и 5 0г , отлична от нуля и мо- 
лекулярное поле имеет в плоскости хі направле- 
ние, в общем случае не совпадающее с направлени- 
ем осей (11-13). Используя ту же процедуру, полу- 
чим нелинейные уравнения движения, соответ- 
ствующие волновому вектору д 

т88 х =и о 8^+І 0 (8 0х ) 2 }88 у + 

+{(/- + 2І-8°*)881 + <Ю- }881 , (26) 

ісо88І = {щ-Ѵ- +2І-8°*)8 0х }88- 2 - 
-{Л5°-'+/ 0 (5°-') 2 }5^- 

-8°*Х1- -{(У ? + 2І-8°*)88; + 8П-}882, (27) 

ісо881 = -со 0 88і. (28) 

Однородная система линеаризованных уравне- 
ний (26-28) при сЮ ? = 0 имеет нетривиальное реше- 
ние только в случае, если определитель этой систе- 
мы тождественно равен нулю. Решая секулярное 
уравнение аналогичным образом, получим в пред- 
переходном состоянии 

«,(?) = 0, (29) 

а два других имеют вид 

соІ(д) = Ѵ 0 8 0: +І,(8 0: ) 2 ] 2 + 
+со 0 (со 0 -^8 0х -2І-8 0х 8 0х ). (30) 

Поскольку для исходной (конечной) структуры 
при Т— > Т е выполняется 5 ,,к — »со 0 // 0 , то моды псевдос- 
пиновых волн со ( = 0 , а а> 2 , з является критической 
только в длинноволновом пределе ( 7 — >0 (29, 30). В 
этом предельном случае квадрат частоты мягкой мо- 
ды определяется спонтанной разностью заполнения 

® 2 2 , (?) = Ѵо8°' + 1ц (^ ) 2 ] 2 - ( 2 /, / ^)со;8* .( 31 ) 

При Г— »Г С и Т С -^Т Ш параметр порядка 5 % стре- 
мится к нулю [3], следовательно, собственные час- 
тоты со 23 (0) также стремятся к нулю (31). В то время 
как в предпереходном состоянии мода псевдоспи- 
новых волн может стать мягкой модой при любом 
значении критического волнового вектора (из об- 
ласти 0 << 7 0 ^? ад г), для которого имеет максимум, в 
исходной (конечной) структуре мода становится 
мягкой только при ( 7 — >0. Если критическое значе- 
ние волнового вектора в предпереходном состоянии 
равно 77 -> 0 , тогда собственные векторы флуктуаций 
псевдоспинов одинаковы как в предпереходном 
состоянии, так и в исходной (конечной) фазе. 


4. Уравнения динамики системы псевдоспинов 
с учетом релаксации 

Перепишем для удобства векторное ур. (6) в виде 
системы трех уравнений для компонент псевдоспина 8 

8 х = [/ 0 5- _ + 1 0 (8 : ) 2 + П]8 } - у[5“ - 5“' ], (32) 

8 У = -\.І„8 : + / 0 ( 5 : ) 2 + 0]5 - со 0 8 х - у8 ѵ , ( 33 ) 

8 : = со 0 8 у -у{8 : + Ѵ 0 8 Х + 1 0 (8 : ) 2 + Щ8 * со 0 '}. (34) 

Покажем, что из (32-34) в низкочастотном ре- 
жиме вытекает уравнение типа Ландау-Халатнико- 
ва. Дифференцируя по времени и используя прос- 
тые алгебраические преобразования, сведем 
(32-34) к системе двух уравнений: 

5-' + у[ 2 + со~'8 0х (У 0 + 2І п 8 І1г )]5* + 

+ую 0 " 1 5 0д О + К 2 + у 2 ](1 + а- '8 0х ^ 0 )8 :: + 
+а>- 0 і 8 0х І 0 8 0х [<о 2 о+Г 2 ](81 2 + 

+со 0 Ѵ 0 8 г + І 0 (8*) 2 +П]88 х =-(о 0 8 а: [б) 2 0 + у 2 Г1, (35) 

88 х = -у88 х + ю 0 1 [7 0 5 г + 1 0 (8 : ) 2 + О] х 
х(8 х +у{8 : +8 0х а 0 Ѵ 0 8 г +І 0 (8 : ) 2 +П]}), (36) 

где 88 Х =8 Х -8° Х . Разрешая систему ур. (35, 36) отно- 
сительно О как функции параметра порядка и 
его производных методом последовательных приб- 
лижений по дисперсии, диссипации и нелиней- 
ности, получим во втором порядке: 

Г + уГ +&]8? -(ѴЛЖ 2 + Г 2 ]х 

х(8 х ) 2 +(1/2)(8 х ) 3 = ё 2 сХ П, (37) 

х = х 4 /(і- ^ х е/ ) = (Л(1 + )}■', 

® 2 =К 2 + у 2 ](1 + й)- 1 / 0 5°Д. (38) 

С другой стороны, согласно уравнению Ландау- 
Халатникова [5], можно записать 

8 х + Г 8 г =-ѵ8Р /88 х , (39) 

где ѵ - коэффициент, Р - свободная энергия систе- 
мы псевдоспинов в механическом поле, которая 
вблизи Т с имеет вид [3] 

Р = Р 0 +кІ 0 {( 1 / 2)[1 - ф(П(й 0 / 2 к в Т)] х 
х(5Д 2 -(1/3)(/ 0 // 0 )(5Д 3 + 

+(1 / 4 ){Г- / 2со 2 )[1 - ц](8 х ) 4 } - Ш8 Х , (40) 

Л = {\I2){^-^~')\а{^ + \I <7 — 1 ), 

+(1 / 4)(8 2 / 2ю 2 )[1 -і]](8 : ) 4 } - Ш8 Х . (41) 

Из (37-41) видно, что статическая восприимчи- 
вость х имеет полюс при Т=Т„ как того требуют 
общетеоретические положения. Подставляя (40) в 
(39) и приравнивая друг другу правые части (37) и 
(39) получим 

у 2 = й ) 2 [77 /(1 -77)],Я=4ф7/(1-77)] 1/2 , 

ѵ = [со 2 / Й/ 0 (1-77 )], (42) 
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Естественные науки 


® 2 = ®0 2 [1-^] _1 (1 + ®0 Ѵх <5 = (43) 

Таким образом, постулировав выполнение за- 
кона сохранения (5), и при наличии релаксации, 
нам удалось получить выражение (42) для диссипа- 
тивного параметра у в виде функций от ю 0 ,/ 0 . 

Как следует из (42, 43), условие незаторможен- 
ности мягкой моды может быть записано в виде 

г) < [1 - дЩТі со 0 / 2 к в Т))]. (44) 

Анализ (44) показывает, что данное условие вы- 
полняется лишь в узком интервале значений пара- 
метра дюі. Отсюда приходим к выводу о том, что в 
окрестности структурного перехода исходная 
структура - предпереходное состояние - конечная 
структура, где ф>>1 (слабое, среднее туннелирова- 
ние [3]), мягкая мода переторможена и динамика 
параметра порядка становится чисто релаксацион- 
ной. Общее исследование низкочастотной динами- 
ки параметра порядка в окрестности структурного 
перехода может быть проведено на основе ур. (37). 
Однако соответствующий анализ является доста- 
точно громоздким. Случай же ф>>1 представляется 
более простым и физически наглядным. В пределе 
#»1, ур. (37), при 0=0 принимает вид уравнения 
Колмогорова-Петровского-Пискунова в простран- 
ственно-одномерном случае: 

=-у[1-ф(Псо 0 /2к в Т)]5? + 

+у (7 0 / / 0 )(5-') 2 - у - 1 (Л 2 / 2)(8Г + °д 2 ?/д2 2 , (45) 

где 1)хуа 2 т(о () а 2 а10 а с 1 (ІО -8 см=10~ 3 см 2 -с~‘ - эффек- 
тивный коэффициент диффузии, а - межатомное 
расстояние в кристаллической решетке, ѵ звук = 
=(й) 0 7)) 1/2 »1О 5 см-с ч -скорость звука в металлах. Как 
отмечалось выше, в низкочастотном длинноволно- 
вом приближении слагаемое, ответственное за 
пространственную дисперсию, может быть учтено 
аддитивным образом, простым добавлением в (37). 

Полученные результаты для перехода исходная 
фаза - предпереходное состояние - конечная фаза, 
вызванного изменением температуры, могут быть 
пролонгированы на случай структурного перехода, 
стимулированного изменением внешней силы. 
Согласно результатам [3] в окрестности структур- 
ного перехода, вынужденного изменением внеш- 
ней силы, поведение параметра порядка имеет ре- 
лаксационный характер, что приводит к связи па- 
раметра порядка с компонентой механического по- 
ля напряжений 0,=7гЯ _1 (<т)^. Поэтому в окрест- 
ности структурного перехода, стимулированного 
деформацией, учет внешней силы сводится к пере- 
нормировке константы двухчастичного взаимодей- 
ствия / 0 =/ 0 +7гЛ _1 (сг), а динамика параметра порядка 
подчиняется ур. (45). 
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5. Обсуждение результатов и выводы 

Одним из главных свойств рассматриваемого 
структурного перехода исходная структура - предпе- 
реходное состояние - конечная структура является 
то, что спонтанное появление отличного от нуля па- 
раметра порядка в исходной (конечной) фазе нару- 
шает внутреннюю симметрию системы. Частота мо- 
ды, восстанавливающей симметрию, должна непре- 
рывно уменьшаться при Т-^Т + (Т ). В случае струк- 
турного перехода почти второго рода при Т=Т + (Т ) 
частота этой мягкой моды практически обращается в 
нуль. С другой стороны, при подходе к Т~( Т ) со сто- 
роны предпереходного состояния частота мягкой 
моды понижающей симметрию предпереходного 
состояния также уменьшается при Г— >Г + (Г~). 

Восстанавливающая сила для смещений атомов, 
отвечающих такой моде, стремится к нулю до тех 
пор, пока псевдоспиновая волна не конденсируется 
на границе устойчивости. Следовательно, статичес- 
кие смещения атомов при почти непрерывном пере- 
ходе из предпереходного состояния в исходную (ко- 
нечную) структуру представляют собой заморожен- 
ные смещения мягкой псевдоспиновой волны. Па- 
раметром порядка при таком переходе является ста- 
тическая компонента собственного вектора мягкой 
псевдоспиновой волны. Так как исходная (конечная) 
структура характеризуется макроскопической спон- 
танной положительной (отрицательной) разницей 
заполнения левого и правого минимума одночастич- 
ного потенциального рельефа, одинаковой на каж- 
дом узле, то мягкая псевдоспиновая волна является 
длинноволновой (1? — >0). Следовательно, псевдоспи- 
новая волна является мягкой модой, с одной сторо- 
ны, понижающей симметрию предпереходного сос- 
тояния, а, с другой стороны, восстанавливающей 
симметрию исходной (конечной) структуры. 

Феноменологический учет релаксации псевдос- 
пинов в окрестности структурного перехода, осно- 
ванный на выполнении закона сохранения квадрата 
длины вектора Блоха, который подчеркивает коллек- 
тивный, когерентный характер релаксационных про- 
цессов, приводит к следующим результатам. Во-пер- 
вых, удалось подтвердить такие известные эффекты, 
как смягчение коллективной моды, критическое за- 
медление. Во-вторых, из предложенного подхода ес- 
тественным образом вытекает важный вывод о том, 
что мягкая мода в окрестности структурного перехо- 
да исходная структура -предпереходное состояние - 
конечная структура переторможена и, следовательно, 
динамика параметра порядка имеет релаксационный 
характер. В-третьих, поскольку релаксационные про- 
цессы приводят к связи компоненты механического 
поля напряжений с параметром порядка [3], дина- 
мика параметра порядка в деформируемом кристалле 
также имеет релаксационный характер. 
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